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Abstract 

Recently, is has been observed that events with the same total transverse energy of light charged 
particles (LCP) in the quasi target region, Ej_^2-i show two quite distinct reaction scenarios in 
the projectile domain: multifragmentation and residue production. This phenomenon has been 
dubbed " bimodality" . Using Quantum Molecular Dynamics calculations we demonstrate that this 
observation is very general. It appears in collisions of all symmetric systems larger than Ca and at 
beam energies between 50 A.MeV and 600 A.MeV and is due to large fluctuations of the impact 
parameter for a given Ej_^2- Investigating in detail the Ej_^2^™' i'^ which both scenarios are present, 
we find that neither the average fragment momenta nor the average transverse and longitudinal 
energies of fragments show the behavior expected from a system in statistical equilibrium, in 
experiment as well as in QMD simulations. On the contrary, the experimental as well as the 
theoretical results point towards a fast process. This observation questions the conjecture that the 
observed bimodality is due to the coexistence of 2 phases at a given temperature in finite systems. 

PACS numbers: 24.10.Lx, 24.60.Lz, 25.70.Pq 
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I. INTRODUCTION 

A while ago the INDRA and ALADIN collaboration has discovered [l| that, in collisions 
of heavy ions - Xe+Sn and Au+Au between 60 and lOOA.MeV incident energy in a 
small interval of the total transverse energy of light charged particles {Z < 2) in the quasi- 
target (QT) domain, -E^]^, a quantity which is usually considered as a good measure of 
the centrality of the reaction, two distinct reaction scenarios exist. In this Ej^i2^nteTYal, in 
forward, quasi-projectile direction, either a heavy residue is formed which emits essentially 
light charged particles, or the system fragments into several intermediate mass fragments. 
In the original publication this phenomenon has been termed "bimodality" due to reasons 
we will discuss further on. We stick to this name although our interpretation of the origin 
of this phenomenon is different, and we will call this interval in i^j^^n" bimodality interval". 
In the meantime, this effect has also been observed by other groups 2|. 

This observation has created a lot of attention, because a couple of years before theory 
has predicted j^, 4, 5] that in finite size systems a first order phase transition weakens: in 
a finite size canonical ensemble, which is determined by the temperature T, the number 
of particles N and a given volume V (or a given pressure p), it becomes more like a cross 
over. In infinite matter the two phases coexist only at the transition temperature. Below 
the transition temperature, Tt{N, V), the system is in one phase and above Tf in the other 
phase. Because energy fluctuations are suppressed by oc this statement is also true 

if the large system is treated micro-canonically. In finite systems, the situation is different. 
In a canonical description of the system, for a given T, N and V, the energy fluctuations can 
become large, even larger than the finite latent heat. Therefore - for a given temperature 
close to Tt - for the same values of T, N and V, the system can either be in the gas or 
in the liquid phase. This means that if the system stays for long in thermal equilibrium, 
it moves back and forth from one phase to the other. The simultaneous appearance of 
these two modes (phases) for a given value of T,N,V has been called "bimodality". In a 
micro-canonical description, bimodality is not possible, because the energy of the two phases 
differs. 

It is of course all but easy to identify these theoretical results with observables obtained 
in a heavy ion reaction. Assuming that -E^^is ^Iso a measure of the temperature of the 
system it is nevertheless tempting to identify the residue with the liquid phase of nuclear 
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matter, and the creation of several medium or small size fragments with the gas phase. The 
experimental observation of the above mentioned bimodality scenario would then just be 
the experimental confirmation of the theoretically predicted bimodality. 

If this were true, the observation of bimodality would solve a longstanding problem of 
heavy ion physics, the quest to identify the reaction mechanism which leads to multifrag- 
mentation. This problem arrived because many observables could be equally well described 
in thermodynamical or statistical theories [y, |7| as in molecular dynamics type models 



lid, although the underlying reaction mechanism is quite different. The statis- 
tical models assume that the system is in statistical equilibrium when its density reaches 
a fraction of normal nuclear matter density. Then, it suddenly freezes out and the frag- 
ment distribution is determined by phase space at freeze-out. In dynamical models, on the 
contrary, fragments are surviving initial state correlations which have not been destroyed 
by hard nucleon-nucleon collisions during the reaction, and equilibrium is not established 
during the reaction. They can be already identified very early in the reaction, when the 
density is still close to nuclear matter density. A detailed discussion of how the reaction 
proceeds in these models can be found in 

Recently, it has been shown that the observation of bimodality alone does not allow for 
the identification of the reaction mechanism. Dynamical models describe the bimodality 
signal as well. In the bimodality i^^^interval, they also show the presence of two different 
event classes, and reproduce quantitatively the scaling properties of -E^^ 13]. Therefore, 
further studies of the bimodality interval in Ej_^2^ie necessary to elucidate the reaction 
mechanism. 

The variables T, N and V determine the canonical ensemble completely. If bimodality, 
in the sense of the coexistence of the two phases at a given temperature, is at the origin of 
the experimental observation, the values the average source velocity of both modes as well 
as their temperature have to be identical. 

In this article, we analyze the experimental and theoretical events which fall in the bi- 

OT' 

modality Ej_]^2^Taiei\aX, without further cuts, to study the average system properties in this 
interval, and to investigate why bimodality, defined as above as the observation of two dif- 
ferent reaction scenarios in a narrow interval of -E^]^, is feasible even if the system has not 
reached thermal equilibrium. In section II, we investigate in detail the experimental ob- 
servables in the bimodality £^^^bin and whether they are compatible with the assumption 
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that there are two phases in thermal equihbrium. This detailed study is possible due to the 
very good acceptance properties of the INDRA detector. In a second step, we investigate in 
chapter III whether bimodality is a phenomenon which occurs only in a small range of sys- 
tem sizes and beam energies in which the system hits the transition temperature, or whether 
it is a more general phenomenon. Especially, the energy dependence is of interest, because 
with increasing energy there is a change in the type of matter from which fragments are 
formed. At low energies, it is the participant matter (the overlapping part of projectile and 



target) which forms the fragments, whereas at energies of a couple of hundreds o 



the fragments are formed from spectator matter (the non overlapping part) 12, llSj- 



A.MeV, 



II. SYSTEM PROPERTIES IN THE BIMODALITY INTERVAL OF E' 



QT 
±12 



For the investigation of the physics in the bimodality E'^^^interval, we follow the defini- 
tions of ref. 1]. -Ej^i^is defined as the total transverse energy of particles with charge Z < 2 
on the quasi target side {6cm > 90°), calculated in the frame in which the momentum tensor 
of all fragments with charge Z > 3 is diagonal. The diagonalization is done event by event. 
We define 02 as 

a2 = (Zmaxl ~ Zmax2) / (Zmaxl + "Z'max'l) (1) 

where Zmaxx is the charge of the largest fragment, while Zmax2 is the charge of the second 
largest fragment, both observed in the same event in the quasi-projectile (QP) hemisphere 
- at polar angles O^m < 90° - in the center of mass of the system. For a more accurate 
extraction of 02, we reject events where less than 70% of the charge of the projectile has 
been detected. Bimodality means that there exists a narrow interval in E'^^in which events 
with large and small 02 values are observed. In this narrow transition region, we expect 
two types of events: One with one big projectile residue accompanied with some very light 
fragments (large 02), the other with two or more similarly sized fragments (small 02). Events 
with intermediate values of 02 should be rare. In Il3[ we have studied function of 

i^^^^for the system Au+Au between 60 A.MeV and 150 A.MeV incident energy that has 



been measured by the INDRA - ALADIN collaboration at GSI [ij]. We focus here on the 
data at 60 A.MeV and concentrate on that experimental bimodality interval of E'^^^ where 
the transition from small to large a2 values occurs. We compare the experimental events 
in this interval with filtered numerical simulations for the same Ej^'^^2^aAues. The filtering is 
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done using a software replica of the INDRA experimental set-up. 

The simulations are performed with one of the dynamical models which has frequently 
been used to interpret the multifragmentation observables, the Quantum Molecular Dynam- 



ics (QMD) approach 8|, |9|, ll2|. This approach simulates the entire heavy ion reaction, from 



the initial approach of projectile and target up to the final state, composed of fragments 
and single nucleons. Here, nucleons interact by mutual density dependent two body inter- 
actions and by collisions. The two body interaction is a parametrisation of the Bruckner 
G-Matrix supplemented by an effective Coulomb interaction. For this work, we have used 
a soft equation of state. The initial positions and momenta of the nucleons are randomly 
chosen and respect the measured rms radius of the nuclei as well as the Fermi distribution 
in momentum space. Collisions take place if two nucleons come closer than r = y^cr/vr, 
where a is the energy dependent cross section for the corresponding channel (pp or pn). 
The scattering angle . chosen .andom.y. respecfng the ex_ta..y measured 
Collisions may be Pauli blocked. For details we refer to ref.[8|, 19||. For the later discussion, 
it is important to note that, even for a given impact parameter, two simulations are not 
identical, because the initial positions and momenta of the nucleons as well as the scattering 
angles are randomly chosen. It has been shown [l^ that these simulations give a bimodality 
signal in the same E'^^bin as the experiment and that the beam energy dependence of this 
transition is reproduced. 



A. ai dependence of observables 

If the experimental signal, i.e. the existence of two distinct reaction scenarios in a certain 
E'^^interval, reflects the coexistence of two phases of a thermally equilibrated system, the 
behavior of several observables can be predicted. In this case, both phases must have the 
same source velocity in the cm. and the same temperature. Considering the nuclei resulting 
from the decay of the projectile spectator as an ideal, classical, non interacting gas, the mean 
kinetic energy of each fragment or nucleon in the center-of-mass frame of this gas should be 
3A;T/2. This has also to be true for the largest fragment, called projectile remnant (PR), 
in contradistinction to other theoretical approaches in which the PR properties reflect the 
violence of the reaction, which depends on centrality and hence on its size. The longitudinal 
velocity loss and the transverse velocity of the heaviest fragment have been systematically 
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studied for different projectile-target combinations 16|. The value of both differs from the 
expectation for a heat bath particle, and depends on the number of nucleons the projectile 
has lost in the course of the interaction. This is one of the interests to study first the 
properties of Z^axi in fhe bimodality E'^^'^interval. The other is that observables of the 
heaviest fragment should be give a clear signal of the process because they are the least 
spoilt by eventual preequilibrium processes which may disturb the light particle spectra. 
Fig. [T]shows the transverse velocity i3± = v±/c (top) and the longitudinal velocity, i3\\ = v\\/c 
(middle), in the cm. of the system, for the heaviest fragment in the QP region as a function 
a2 and for those fragments which belong to events which fall into the bimodality interval, 
12 < £^^2/ -^ocm < 17. Eocm IS the incident energy per nucleon in the cm. system. The 
bottom panel displays the average charge of the heaviest fragment as a function of 02- In 
all figures, we show the INDRA data as points and compare the experimental results with 
unfiltered (full line) and filtered (dashed line) QMD predictions. We see that the filter 
changes substantially /3± (due to a cutoff at small transverse velocities) and to a less extend 
l3\\. After filtering, the QMD calculations reproduce the trend of the data but overpredict 
the transparency, i.e. underpredicts (3±, for small 02 events corresponding to those with a 
small Zmaxi- If our assumption that E'j^^gch^r^cterizes the excitation energy and hence the 
temperature of the system and that the events in the bimodality bin have that temperature 
at which the system can either in the gaz or in the liquid phase we expect that l3± as well as 
P\\ -Psource oc l/\/mass. Here Psource is the velocity of the bimodal system in the cm. which 
should be independent of 02 because the velocity of the bimodal system does not depend on 
which mode is realized. 

Consequently, this assumption does not offer the possibility that (3\\ increases with the 
particle mass or charge. The experimentally observed linear increase of l3\\ with the mass 
of the fragment follows, however, the trend already observed in central and semi-central 
Xe+Sn events at 100 A.MeV [l^ (where no bimodality has been observed). The average 
velocity of the heaviest fragment, in events with a smaller E'j^^gthan at bimodality follows 
this systematics as well. This linear decrease of the PR velocity with decreasing mass is a 

for higher beam energies by Ricciardi et al. [17|, although it has not been shown yet that this 
systematics is still valid for such small PR in this energy domain. Such a linear dependence 
is expected if nucleons are removed randomly from the cold projectile nucleus, under the 
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FIG. 1: (Color online) Average transverse velocity = u_l/c (top) and /3|| = fy/c (middle) of the 
QP fragment {dcm < 90°) with the largest charge, Zmaxi, as well as Zmaxi (bottom) as a function 
of the QP 02 for those events which fall in the bimodality interval -£'_li2' text. We compare 
INDRA data (symbols) with filtered (dashed lines) and unfiltered (full lines) QMD simulations for 
Au+Au collisions at 60 A.MeV bombarding energy. 
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condition that they do not interact with the residue anymore. This removal leads to a 
deceleration and an excitation of the remnant. Thus this increase of can be understood 
in models in which the heaviest fragment is not in thermal equilibrium with the emitted 
particles, but finds no explanation in purely thermal models. 

In addition, we observe that the average transverse velocity decreases with In order 
to discuss the compatibility with the assumption of a canonical bimodal system, we have 
to consider the unfiltered QMD events. They show an almost constant transverse kinetic 
energy of 35 MeV, independent of the fragment mass and charge, as it will be discussed at 



the end of section II. B. This value confirms the analysis of [18| and is too large (even if one 
considers a radial fiow, as we discuss in the next section) to be compatible with kT: the 
expected kinetic energy for fragments in a thermal heat bath has to be smaller than the 
binding energy of nucleons. Hence, in the bimodality i^^^interval, also is incompatible 
with the assumptions that two phases are present, having both the same temperature. 

Thus neither the mass dependence of P\\ nor that of P\\ of the heaviest fragment Zmaxi 
in the bimodality Ej_^2^ntei\a\ are compatible with the expectations for a finite system in 
which two phases are in equilibrium. On the contrary, they follow the systematics which 
we have observed for other E'];^ 2 regions where no experimental signs of the presence of two 
phases in thermal equilibrium, i.e. bimodality behavior, are found. These properties can be 
explained in models which are genuinely non-equilibrium and which have successfully been 
applied to interpret data in many experimental situations. 

B. Particle properties in the E'j^^gi'^terval which shows bimodaUty 

A thermal system has to be isotropic in its rest system. The system we study here is the 
ensemble of all QP particles, i.e. those which are emitted at angles 9cm < QC. The degree of 
isotropy can be studied with help of the momentum tensor in the rest system of the source 

_ 2<pi>-<^>-<pi> 

Pz is the momentum in the beam direction. Qzz = if in the rest system of the source 
the distribution is isotropic. Negative values indicate a preferred emission in transverse 
direction. Fig|2] gives an overview over Qzz of light charged particles (Z=l,2) and fragments 
{Z > 3) in the bimodality E^^^interval. 
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FIG. 2: (Color online) Distribution of Q^z for Au+Au events at 60 A.MeV which fall in the 
ii^_£j^2biniodality bin. We display the Qzz distribution (see text) for light charged particles (Z=l,2) 
(top row) in their own center-of-mass, and for fragments with charge Z > 3 in their own (middle 
row) and in the Z=l,2 cm. system (bottom row), as a function of 02. The points are the mean 
values of Qzz as a function of 02. ^From left to right, we display unfiltered and filtered QMD 
predictions, and the INDRA data, respectively. The points mark the mean values. 
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FIG. 3: \Qzzi foi^ Au+Au at 60 A.MeV in the £'_£^r,bimodality bin, calculated for the second largest 
charges Zmax2 of all events taken all together (around their average momentum), as a function of 
a2. We compare INDRA data (symbols) with filtered (dashed lines) and unfiltered (full lines) QMD 
predictions for Au+Au collisions at 60 A.MeV. 

The average longitudinal velocity of all quasi projectile nucleons in the cm. depends on 
a2. In the bimodality scenario this should not be the case. It is therefore not meaningful to 
analyze all events of the bimodality E'j^j^gi^terval in a common rest frame. We concentrate 
here on the question whether for a given 02 the events are isotropic. This a necessary but - 
as just mentioned - not a sufficient condition for bimodality. To obtain the Q^z distribution, 
we adopt the following procedure: We select events which fall into the iJ^j^interval which 
shows bimodality and calculate for each 02 bin the average velocity of a) all LCP (Z=l,2) 
and b) of all fragments Z > 3 of all events which fall in this bin. In one of the rest systems we 
calculate subsequently for each event Qzz of either the fragments or the LCP. This procedure 
is applied to eliminate the dependence of the velocity of the rest system of the fragments on 
02 (see fig. [1]), because we expect isotropy only in the rest system. /^From top to bottom, 
we display the Qzz distribution of the LCP in the rest system of the LCP, that of fragments 
in the rest system of the fragments and that of fragments in the rest system of the LCP. 
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From left to right we show the unfiltered, the filtered QMD predictions, and the INDRA 
data, respectively. Please note that the Ej_T^2^\iieiYa\ is determined by the energy of the 
Z=l,2 particles in the QT domain, whereas here we study the properties of Z=l,2 particles 
in the QP domain. Thus, autocorrelations between i^^^^and the particles studied in the QP 
domain are minimized. 

We see from the top figures that the LCP in this E'^^interval are preferable emitted into 
transverse direction in their rest frame {Qzz < 0). The anisotropy depends only slightly on 
02- It may be due to the fact that not only the bimodal system but also preequilibrium 
emission contributes to the spectra. It indicates that the origin of the Z=l,2 is never a pure 
thermal source, neither for small 02 nor for large 02 events. The experimental filter changes 
little as far as the LCP anisotropy is concerned and the filtered QMD predictions agree quite 
well with the INDRA data. In their proper rest system the fragments are emitted preferably 
in forward/backward direction. The experimental filter brings the QMD predictions of the 
average Qzz value closer to zero, in agrement with the experiment. At very large values of 
02, the emission becomes isotropic. Such a isotropy is expected for example for the emission 
of a light charged fragment from a compound nucleus. For small 02 values, Qzz of the 
fragments in the rest system of the LCP fluctuates around as expected for an isotropic 
source. Thus small 02 events come closest to isotropy. For 02 > 0.5, Qzz is close to two. 
This means that in the LCP rest system, the largest fragment is preferably emitted in the 
beam direction. This observation rules out, on the other hand, the hypothesis that the 
ensemble of fragments and light charged particles at large 02 values can be considered as 
a pure liquid phase. This would require that the direction of the velocity of the largest 
fragment is randomly distributed in the rest system of the liquid. 

One may argue that the largest fragment is not really in thermal equilibrium in the 
sense that it has not lost completely its memory on the entrance channel, i.e. its initial 
velocity direction. Such an argument cannot be put forward for the second largest fragment. 
Therefore, we display in figj3]the average Qzz for the second largest fragments calculated in 
the rest system of those fragments. The INDRA data (symbols) as well as the unfiltered 
(full line) and filtered (dashed line) QMD events show Qzz < independent of 02. The 
value of Qzz ~ —0.2 for large 02 indicates that the emission of the second largest fragment 
is essentially random but still preferably in transverse direction in the fragment rest system. 
For the small 02 values we find Qzz ~ —0.4. Again, the second largest fragment is preferably 
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emitted in transverse direction. For the interpretation, we have to combine this result with 
that shown in the lower right panel of fig. [2j in the rest system of the fragment, the 
emission of all fragments is almost isotropic. A small value of 02 means that the two biggest 
fragments have about the same size. Therefore, combining Q^z of all fragments and of the 
second largest fragments yields the following scenario: For small 02 values the fragment 
with the highest charge has the largest velocity in beam direction, but this time the lighter 
fragment can balance the momentum and therefore the total emission pattern appears to 
be isotropic, although, if one looks into the detail, it is not. 

Thus, only for the largest 02 values, the fragments are isotropically distributed in their 
rest system as they should if the QP system represents the vapor or the liquid phase of a 
system in thermal equilibrium. This observation for large a2 values is also compatible with 
the emission of a light fragment from a compound nucleus. 

Fig. m shows the average velocity in longitudinal f3\\ cm and transverse direction /3± in 
the reaction cm. system as a function of the fragment mass, for those events which fall in 
the bimodality interval of -E^^(left), and for all events (right). In all figures we display the 
INDRA data as points and compare them to the unfiltered (full line) and filtered (dashed 
line) QMD predictions. Similar to fig. [1], we observe also here an increase of l3\\ with 
increasing fragment mass. For large fragments the velocity approaches that of the beam 
(/3|l = .179). The fragments in the selected -E^^bin show - as Z^axi - an almost identical 
behavior to those observed without a selection in -E'_li2 5 and follow the Morrissey systematics 
IgI ]. QMD simulations reproduce rather well. The transverse velocity of the heavy 
fragments is small and increases with decreasing fragment mass. Such a dependence has 
been observed for systems with a collective radial velocity 28| . In such a scenario, the form 
of the velocity dependence of the mass allows for the determination of the temperature and 
of the radial velocity if the emitting source is thermal with a collective radial velocity 28|. 
In a thermal system, E± = kT. The dotted curve, calculated according to the formula 
given in ref. 281], gives best agreement with the experimental data for E± = 13.5MeV. This 
value is too high to be compatible with the assumption of a thermalized system. At such a 
temperature, the fragments would not exist anymore. 

On the contrary, the observed value of 13.5 MeV is well described by the Goldhaber model 
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22| . It is based on the assumption that nucleons or fragments are removed rapidly from 



a cold nucleus and therefore their momentum distribution is given by the Fermi motion. 
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FIG. 4: (Color online) Average longitudinal < > and average transverse velocity < /3_|_ > in 
the reaction center of mass system as a function of the fragment mass. In the left panels, we 
display these quantities for those fragments which belongs to events in the bimodality interval of 
-^±12' ^^'^ right panels for those of all events without any iiJ^^cut. INDRA data (symbols) 

are compared with filtered (dashed lines) and unfiltered (solid lines) QMD predictions for Au+Au 
collisions at 60 A.MeV. The dashed bold line is the result of the fit of the experimental data, using 
the approach of ref. 



Then, the average transverse squared momentum per nucleon of a fragment of size Ap is 
given by the Goldhaber formula 22l |: 



pAA.) - 3 ■ -,PfA.^^ (3) 
where Ap[Ap) is the projectile (fragment) mass and pp denotes the Fermi momentum. 
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Therefore 



E^[Af) = 4^ ^ ^ UMeV (4) 



is almost independent of the fragment size and in agreement with the data. The QMD simu- 
lations reproduce the form and the absolute value of (3±{A). This is not astonishing, because 
in this model, fragments are surviving initial correlations which preserve approximately the 
transverse momentum they had initially, and therefore, the average transverse energies of 
intermediate mass fragments are, independent of the impact parameter, close to the value 
expected from the Goldhaber formula. Thus, the fragment average transverse velocities are 
understandable if one assumes that there is a collective radial expansion of the system which 
is superimposed to an average transverse energy, given by the Goldhaber model. 

Fig. [5] displays the experimental average transverse energy, as a function of the 
fragment mass of QP products. We see indeed that E_[_ in unfiltered QMD predictions is 
almost independent of the fragment size. Filtering for the experimental acceptance increases 
E_[_ for intermediate mass fragments and brings the calculation closer to the experimental 
observation. E\\, on the contrary, is strongly mass dependent in the rest system of the QP. 
Up to A=50, this dependence is well reproduced by QMD predictions. Above, there are 
discrepancies. 

A thermal system has to be isotropic in coordinate and momentum space, and for each 
degree of freedom, the average kinetic energy of the fragments has io he E = ^kT in the 
rest system of the source. Therefore it is meaningful to calculate the deviation of 

2<(/3||-<^ll>P> 

from one. The INDRA data as well as the QMD predictions, displayed in fig. [5l show that R 
strongly depends on the fragment size. Up to mass A=40 , R decreases strongly and increases 
slightly above. Therefore, this specific bimodality Ej^i2^nteTYal shows the same behavior 



which has been observed in 



231] for central collisions, where bimodality does not occur. 



The discrepancy above A=40 between the experimental data and the QMD predictions 
disappears if one rejects the symmetric and asymmetric fission events (by requiring that the 
product of the two largest charges is smaller than 600). Thus, these deviations come from 
events in which two fragments of similar charge are observed. 
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FIG. 5: (Color online) Twice the average longitudinal (top right) and the transverse (top left) 
kinetic energy of QP products in the QP ceilfpr of mass system, as well as the ratio of both 
(bottom), as a function of their mass. We compare INDRA data with filtered (dashed lines) 



III. SYSTEM SIZE DEPENDENCE OF BIMODALITY 



After having established that most of the data for Au+Au colhsions at 60 A.MeV are 
quantitatively described in QMD calculations, we study now the system size dependence of 
the bimodality. Unfortunately, no data have been published so far to verify these theoretical 
predictions. In order to discuss the physics, we present, in fig. [6l 02 as a function of the 
reduced -E_li2 (right) as well as a function of the reduced impact parameter h/h^ (left). Here, 
E±i2 is calculated over all particles with charge Z=l,2 (QT and QP) and 02 is given for 
QP products. To account for system energy and size scalings, -E'_li2 is divided by the total 
mass of the system and by the energy per nucleon in the center of mass. The top row 
shows a2 calculated with all fragments. For the bottom row, we require that the fragment 
with the largest charge, Zmaxi-, and that with the second largest charge Zmax2 have both a 
charge larger than 2. We start the discussion with the Au+Au reaction (top left) for which 
we zoom on the narrow impact parameter interval 5 fm < b < 7 fm, where bimodality 
occurs. If plotted as a function of E±i2 for this impact parameter interval, we clearly see 
the two event classes with a distinct 02 value for the same small Ej_i2 interval. Whether we 
limit the events to Z^axi > 2 and Zmax2 > 2 or not does not make a difference, because 
there are almost always two fragments with Z > 3. When plotted as a function of b/bo, 
the bimodality structure with a sudden jump has disappeared. There, the events with a 
small 02 value are distributed over a broad range of impact parameters 12]. Because the 



fluctuations in E±i2 for a given impact parameter are large, some of these events appear in 
the same bimodality interval in Ej_i2 as the events with a large 02- This is shown in fig 13 
which displays the filtered INDRA Ej_i2 distribution and that of b/bo for a given unfiltered 
E_i_i2 in QMD simulations of Au+Au reactions at 60 A.MeV. For an incident energy of 100 
A.MeV, the distributions are similar. The bars mark the standard deviation. We see that 
there is a strong correlation between these observables. In particular, the INDRA set-up 
provides a mean linear response - no saturation - to E±i2 over a large range of E±i2 (hence 
to the multiplicity and to the energy of particles). The distributions are, however, quite 
broad, and hence, for a given experimentally measured -E_li2 value, the unfiltered -E'_li2 's as 
well as the impact parameters show large fluctuations. 

Like Au+Au, the smaller Xe+Sn system (figj6l top right) exhibits two distinct maxima 
of a2, with a sudden jump of the most probable value of 02 (depicted by the dashed lines). 
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FIG. 6: (Color online) Double differential reaction cross section (linear color scale, normalized 
to the number of events) as a function of 02 and of the reduced impact parameter b/bo and the 
reduced E±i2 (see text), respectively. The symbols show the most probable value of 02- The dotted 
line is to guide the eye. In the top row of each panel, all events are accepted, in the bottom row 
only those with Zmaxi > 2 and Zmax2 > 2 are shown. We display Au+Au (top left), Xe+Sn (top 
right), Kr+Kr (bottom left) and Ca+Ca (bottom right) reactions at 100 A.MeV incident energy. 
For Au+Au, note that a narrow selection in impact parameter (5/m < b < 7fm) has been applied 
around the bimodality region. 
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FIG. 7: (Color online) Double differential reaction cross section (logarithmic color scale in arbitrary 
units), for Au+Au at 60 A.MeV QMD predictions, as a function of E±i2 divided by the kinetic 
energy per nucleon in the center of mass of the reaction system for INDRA filtered events (top) 
and as a function of the reduced impact parameter b/bo (bottom), respectively, and as a function 
of the unfiltered reduced E±i2 ■ The symbols show the mean values of Ej_i2 /Eocm and the error 
bars the rms of their distribution. 

The two event classes are also seen if the events are plotted as a function of b/bo, and 
it is visible that they are both associated with quite different impact parameters. Thus, 
nuclei disintegrate in two quite distinct pattern, but they belong to quite different impact 
parameters, and hence to quite different reaction geometries. For the even smaller Kr+Kr 
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system (figlHl bottom left), there is still a sudden jump of 02 as a function of E±i2 , but the 
two event classes become less distinct, the relative yield of intermediate 02 values getting 
higher. This can be explained by the fact that the absolute value of -£^^12 is reduced and 
hence the relative fluctuations around the mean value increase for a given 6/60 • Therefore, 
we find again the existence of the two maxima in 02 for (almost) the same value of the 
reduced E±i2 if Z^axi > 2 and Zmax2 > 2. The impact parameter of the two event classes 
is, however, quite different, the minimum of a2 is less pronounced. In addition, events with 
Zmaxi > 2 and Zm.ax2 > 2 becomc rare. The majority of events with a large 02 value are 
now those in which the second largest fragment charge has Z = 2. For such small systems, 
already a beam energy of 100 A.MeV makes the reaction that violent that in central collisions 
fragments do hardly survive. Finally, for the very small system Ca+Ca (figEl bottom right), 
bimodality becomes almost impossible, many intermediate values of a2 are highly populated, 
because the system is too small. For the rare events with Z^axi > 2 and Zmax2 > 2, large 
values of 02 are impossible, and therefore we cannot have two distinct maxima anymore. 
The reaction is dominated by a emission from the QP, as can be seen in the top row. 

It is remarkable that, independent of the system size, the sudden jump of the mean 
a2 value occurs around E±i2 /Eocm/i^p + ^t) ~ 0.2, i.e. when the transverse energy of 
light charged particles per nucleon is identical. This scaling is understandable, because 
this quantity measures the energy transfer in the reaction and extends the scaling we have 
observed already for the beam energy dependence of the bimodality for the Au+Au system 

It is an important observation that nature disfavors intermediate a2 (02 ~ 0.5) values. 
Either a big cluster emits small fragments or we observe multifragmentation, events in which 
several intermediate mass fragments are produced. The dominance of these two reaction 
scenarios is independent of the system size. Dynamical models reproduce this observation 
but it is very desirable to know whether also statistical models predict such a suppression 
of intermediate a2 events. This would allow to elucidate whether this suppression is related 
to phase space or to the nucleon-nucleon interaction. 

Independent of the system size, the events with large 02 belong to other impact parameters 
than those with a small 02. But E^i2 and the impact parameter are not well correlated and 
therefore, due to the fluctuations of E±i2 for a given impact parameter, events with very 
different 02 values appear in the same E±i2 interval. For smaller systems, the events with 
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large 02 are those in which an a particle is emitted from the residue. There we see as well 
that events with different do not belong to the same -E_li2 interval. 



IV. ENERGY DEPENDENCE OF BIMODALITY 




FIG. 8: (Color online) of the projectile spectator as a function of (top) and of (middle) 



for Au+Au collisions at 600 A.MeV [33t |. On the left hand side, we show the Aladin experimental 
results, on the right hand side the QMD predictions. The bottom row displays the theoretical 
predictions of 02 as a function of i?_Li2 (right) which are enlarged around intermediate -E^i2 values 
(delimited by dashed lines) in the left panel. In this latter panel, the symbols show the most 
probable value of 02, and the dotted line is to guide the eye. 

If bimodality is a special manifestation of the general feature that, in heavy ion reac- 
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tions, two distinct reaction scenarios with quite different 02 values exist, it is tempting to 
see whether this observation continues to higher energies. Whereas at low energies multi- 
fragmentation happens in central collisions, and therefore fragments are formed from the 



participating nucleons, at higher energies 



121 ]. multifragmentation happens at large impact 



parameters and fragments are formed from spectator nucleons. There are not many experi- 
mental results available which allow to study this question, especially there is no experiment 
in which simultaneously E^u and all fragment charges have been measured in the target 
or projectile spectator. The only experiments which allow to address this question is the 
ALADIN experiment at GSl. In this experiment, Au + Au at 600 A.MeV incident energy 



33( 1 . almost all fragments with Z > 2 of the projectile spectator (PS), i.e. in forward direc- 



tion, have been measured, but no light charged particles. Because Z = 1 particles are not 
measured, -Ej^fgcannot be extracted, and we have to conclude the existence of bimodality 
indirectly. In ref. 3^, it has been shown, however, that the inclusive (impact parameter 
averaged) yield of 02 has maxima at small and large values, separated by a minimum at 
around a2 = 0.5, similar to the observations at energies around 100 A.MeV. 



These data have also been analyzed by a statistical model approach 31[ where it has 



been shown that the experimental mean values and fluctuations are well described once 
the distribution of the system energy, E, has been adapted. It has, however, not been 
demonstrated that this energy distribution corresponds to that which is expected for a 
given temperature of the system, as required by the bimodality assumption. In any case, 
the experimental inverse slope parameters are much larger than those of statistical model 



calculations [3j] . These experimental parameters are compatible with the prediction of the 
Goldhaber model (eq. [3]). 

In a first step, we have to verify that QMD simulations reproduce correctly the pattern 
of 02 as a function of Zhound, the measured charge of all fragments with 2 < Z < 30 (^62) 
or 3 < Z < 30 (Zbs) in the PS. In the ALADIN experiments, it has been shown that Zbound 
is strongly correlated with the energy deposit during the reaction 33|]. In a second step, we 
replace then Zbound by E±i2 (calculated with both target and projectile spectator particles). 
Because Zbound as well as E±i2 are considered as a good mesure for the centrality of the 
reaction, and more precisely of the energy deposit in the spectator, such a replacement is 
meaningful. Fig. [8] shows in the two top rows a2 as a function of the bound charge in the PS 
domain. On the left hand side, we display the results obtained with the ALADIN set-up, on 
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the right hand side the filtered QMD events. Filtering is here not a very important issue, 
because the ALADIN set-up registers the large majority of the fragments in the PS region. 
02 as a function of Zi,2 and is shown in the top and middle panel, respectively. We see 
a quite reasonable agreement between theory and experiment. This allows us to replace in 
the bottom row Zhound by £'±12 • On the right hand side of the bottom row, we see that also 
at 600 A.MeV bimodality can be observed. In a E±i2 interval around 9 GeV, we see two 
event classes: one with a large and one with a small value of a2- The left hand side of the 
bottom row zooms in this -E'_li2 interval. As at beam energies around 100 A.MeV [3], the 
two event classes are separated by a region with 02 ~ 0.5 which contains only a very limited 
number of events. 



V. CONCLUSION 



The appearance of two distinct reaction scenarios, multifragmentation and residue pro- 
duction for the same value of E±i2 is a very genuine phenomenon in heavy ion collisions. It 
exists over almost the whole energy range for which multifragmentation has been observed, 
and it exists in participant fragmentation as well as in spectator fragmentation. The fact 
that physical events show either a small or a large 02 value but almost never an interme- 
diate a2 value is first of all remarkable. It is also very general, and the classification of 
events with the help of a2 is a good way to elucidate this fact. Whether phase space or 
nuclear interactions are at the origin of the lack of events with intermediate 02 values is still 
unknown. 

In the dynamical QMD model, the large impact parameter fluctuations for a given E±i2 
are the reason that events with small and large 02 values appear for the same value of E±i2 . 
The data are in agreement with predictions of models which assume that multifragmentation 
is a fast process. QMD simulations, in which fragments are surviving initial correlations, 
reproduce the large majority of the experimental observations. 

The investigation of the bimodality E±i2 interval shows that the majority of the events 
in this interval has properties which are not compatible with the assumptions that large 
and small 02 events belong to two phases which exist simultaneously in a small interval of 
the temperature measured by E^i2 . We can, however, not exclude that a subset of the 
events shows the properties expected from a statistically equilibrated system. Investigations 
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whether such a subset can be found have been advanced recently 3J]. It would therefore 
be interesting to see whether such a subset may be a sign of true bimodality, i.e. can be 
reproduced in statistical models with an energy fluctuation expected for a system having a 



fixed temperature. Of course generalized statistical ensembles can be defined 



35| assuming 



that in each event E,N,V differ. It remains to be seen how the distribution of E,N,V can be 



assessed and whether such an approach is compatible with data. 
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